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ABSTRACT 

We determined masses for the 7167 DA and 507 DB white dwarf stars classi- 
fied as single and non-magnetic in data release four of the Sloan Digital Sky 
Survey (SDSS). We obtained revised T cr r and logg determinations for the most 
massive stars by fitting the SDSS optical spectra with a synthetic spectra grid 
derived from model atmospheres extending to log g — 10.0. We also calculate 
radii from evolutionary models and create volume-corrected mass distribu- 
tions for our DA and DB samples. The mean mass for the DA stars brighter 
than g=19 and hotter than T off = 12 000 K is (M) DA ^ 0.593 ± 0.016 M Q . 
For the 150 DBs brighter than g=19 and hotter than T cff = 16 000 K, we find 
(■M)ub — 0.711 ± 0.009A4©. It appears the mean mass for DB white dwarf 
stars may be significantly larger than that for DAs. We also report the highest 
mass white dwarf stars ever found, up to 1.33 M.q- 
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1 INTRODUCTION 

White dwarf stars are the end product of evolution of all stars with initial masses up to 
around 9 A4 Q and their distribution contains information about star formation history and 
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subsequent evolution in our Galaxy. As the most common endpoints of stellar evolution, 
white dwarf stars account for around 97% of all evolved stars. Considering there has not yet 
been enough time for any of them to cool down to undetectability, they can also provide 
independent information about the age of the Galaxy. Through an initial-final mass relation 
(IFMR), we can also study mass loss throughout the stellar evolution process. Because 
white dwarf progenitors lose carbon, nitrogen, and oxygen at the top of the asymptotic 
giant branch, they are significant contributors to the chemical evolution of the Galaxy and 
po ssibly an important sou rce of life sustaining chemicals. 



Kleinman et al. 



stars in the Sloan Dig - 



Eisenstein et al. 



(12006|) 



(120041 ) published the spectra of 2551 white dwar; 
ital Sky Survey (SDSS) Data Release 1 (DR1), covering 1360 deg 2 . 
extended the white dwarf spectroscopic identifications to the 4th SDSS Data Release with 
a total of 9316 white dwarf stars reported, more than doubling the number of spectroscopi- 
cally identified stars (IMcCook fc Sionll2003l ). In both works, the authors fit the entire optical 
spectra from 3900A to 6800A to DA and DB grids of synthetic spectra derived from model 
atmospheres calculated by Detlev Koester, up to log g = 9.0. Their fits include SDSS imaging 
photometry and allow for refluxing of the models by a low-order polynomial to incorporate 
effects of unknown reddening and spectrophotometric errors. The SDSS spectra have a mean 
g-band signal-to-noise ratio SNR(g)^ 13 for all DAs, and SNR(g)^ 21 for those brighter 
than g=19. 

This large sample of stars with spectroscopic fits gives us a new opportunity to fully 
explore the white dwarf mass distribution. Understanding the white dwarf mass distribution 
offers insights into mass loss during stellar evolution, the IFMR, and has bearings on close 
binary star evolution. Our report, as well as many previous studies, detect a substantial 
fraction of low mass white dwarf stars that theoretically cannot have evolved as single 
stars, because the age of the Universe is smaller than their presumed lifetimes on the main 
seq uence. 



Kleinman et al. 



(120041 ) notice an increase in mean logg for stars cooler than T e g 



12000 K, but caution the trend might not be real, indicating a pr oblem in the data or fi t 



technique, instead. The trend has p ersisted into t he larger catalog o 



Madej. Nalezvtv. fc Althausl (120041 ) analyzed the 



Kleinman et al. 



Eisenstein et al 



J2006J). 



(120041 ) sample of fits and 



calculated the corresponding SDSS DR1 pure hydrogen atmosphere (DA) white dwarf mass 
distribution. As expected from the log g trend, they found that the mean mass also increased 
below T e g = 12000 K. Their Table 1 presents all previous mean mass determinations, pro- 
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during an average of 0.57 and a most populated peak at 0.562 JH Q for the 1175 stars 
hotter than T e ff = 12000 K. They did not study the potentially h ighest mass stars with 



logg > 9, because they were limited by the stellar atmosphere fit by 



Kleinman et al. 



(12004] ) 



which artificially pegged stars near the upper logg = 9.0 boundary to the boundary itself. 



The incr ease in mean masses fi tted from optical spectra bel ow T P ^ = 12000 K has been 



seen p rior to 



Kleinman et al 



(120041 ) and has been discussed since iBergeron. Wesemael. &: Fontaine 



(Il99ll ) and iKoesterl (Il99ll ). It is usually dismissed as due to problems in the models: either 
convection bringing up subsurface He to the atmosphere, increasing the local pressure, or 
problems with the treatment of the hydrogen level occupation probability. The new larger 
SDSS data set, however, now gives another opportunity to explore this trend and evaluate 
its cause. 



Most reported white dwarf mass determin ations have been deriv ed b y comparing the 



Kleinman et al. 



feOQJ) and 



Eisenstein et al. 



optica l spectra with model atmospheres, as with 
(120061 ) . For the DA stars, the H7, H8 and H9 lines, in the violet, are the most sensitive 
to surface gravity because they are produced by electrons at higher energy levels, those 
most affected by neighboring atoms. However, these lines are also in the region where the 
atmospheric extinction is the largest and typical CCD detectors are the least sensitive. As a 
consequence, most studies used only the line profiles in their fits, avoiding the dependence on 
often uncertain flux calibrations. The SDSS white dwarf spectra have good flux calibration 
and acceptable SNR redwards of 4000A. The published SDSS catalog therefore fits the entire 
optical spectrum, and not just the H lines, as has been traditionally done. The rationale 
for this approach is the good, uniform spectrophotometry and corresponding broad band 
photometry that can be used in the fits. In addition, a low-level re-fluxing is allowed to take 
out large errors in spectrophotometry and any unknown reddening effects. 



In this paper, we will compare t he measured white dwarf mass distributions from 



Kleinman et al. 



(120041 ) and 



Eisenstein et al 



(120061 ) with previous determinations and attempt to assess the 
reason for the observed increase in mass for lower temperatures. We will also explore the 
observed mean masses and analyze the two different fitting techniques: line profile vs. whole 



spectrum, to see the effects on the resulting mass distributions. 
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2 DATA AND MODELS 

The SDSS imaged the high Galactic latitude sky in five passbands: u, g, r, i and z, and 
obtained spec tra from 3800A t o 9200A with a resolution of ~ 1800 using a twin fiber- fed 



spectrograph (jYork et al 



20001 ) . Since we are primarily interested in the mass distribution 



here, we sele cted only the sin 



6.5 from the 



Eisenstein et al. 



e DA and DB stars with log (7 — cri og5 ^ 8.5 and logg + cri og5 ^ 



(120061 ) sample and re fit them with an expand ed grid of models 



see below) , using the same auto fit routine as in 



Eisenstein et al. 



described in 



Kleinman et al. 



(120061 ) and thoroughly 



(120041 ) . We excluded all stars classified by 



Eisenstein et al. 



( 120061 ) as having either a detectable magnetic field or a companion, metal lines, DABs, and 



DBAs. 



used by 



Our model grid (|Finlev. Koester. fc Basri 



Eisenstein et al. 



1997 



Koester et al. 



200 ll ) is similar to that 



( 120061 ). but extended in T eS and \ogg (100 000 K < T eff < 6000 K, 



10.0 ^ log g ^ 5.0) a nd denser. We chose t he ML2/a = 0.6 parameterization for convection 



as demonstrated by 



Bergeron et al 



(Jl995|) to give internal consistency between tempera- 



tures derived in the optical and t he ultraviolet, photometry , parallax, and with gravitational 



redshift. ML2 corresponds to the 



Bohm &: CassineUil (Il97ll ) description of the mixing length 



theory and a = £/Xp is the ratio of the mixing length to the pressure scale height. The 
models include the and H2 quasi- molecular opacities and only Stark (jLemke!ll9971 ) and 



or < 8500 K. 



Liebert. Bergeron, fc Holberg! (120051 ) 



Doppler broade ning, so the line profiles are not pre cise 

Even though lNapiwotzki. Green, fc Safferl (Il999l ) and 
discuss the necessity of using NLTE atmospheres for the stars hotter than 40 000 K, all 
quoted values are from LTE models, as they also show the NLTE corrections are not domi- 
nant, and our number of hot stars is small. 

To calculate the mass of each st ar from the T P ff a nd log q values obtain ed from our fits, 



we used the evolutionary models of IWoodl (jl995l ) and lAlthaus et al.l (120051 ) with C/O cores 
up to log g = 9.0, and O/Ne cores for higher gravities, Mn e = 10 -2 M*, and Mr = 10~ 4 M*, 
or Mh = 0, to estimate stellar radii for DAs and DBs, respectively. The radius is necessary 
to convert surface gravity to mass. 



3 ANALYSIS 

Before exploring the mass distributions, we wanted to examine the different fitting techniques 
used in the available data sets — the traditional line profile technique and the SDSS whole 
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spectrum approach. We therefore simulated spectra with differing SNRs by adding random 
noise to our models and fit them with our own set of both line profile and whole spectrum 
fitting routines. Our Monte Carlo simulations show that in the low SNR regime, SNR^ 50, 
fitting the whole spectra and not just the line profiles gives more accurate atmospheric 
parameters, as long as the flux calibration or interstellar reddening uncertainties do not 
dominate. We estimate an uncertainty of around AT e g- ~ 500 K and A log g ~ 0.10 at 
SNR=40 for the whole spectra fitting. For SNR=20, similar to the average SDSS spectra for 
g < 19, our simulations indicate AT e g ~ 750 K and A log g ~ 0.15. We do not report in this 
paper on the mass distribution for the stars fainter than g=19 because their smaller SNR 
lead to large uncertainties. Our simulations did not indicate systematic trends between the 
two approaches. 
Although 



Kleinman et al. 



(120041 ) and 



Eisenstein et al 



(120061 ) compared their fits' in- 



ternal errors by fitting dupli c ate sp ectra, they did not display their results as a function 



of temperature. 



Kepler et al. 



(120051 ) specifically analyzed 109 duplicate spectra SDSS DAs 



with 13000 K T e fj 10000 K, near the region where the fit loggs start to increase. They 
showed that the mean fit differences were ax cS — 300 K and eri og9 ~ 0.21 dex for the same 
object but different observations. These values a re larger than the inter nal u ncertainty of the 



Kleinman et al 



2004) 



Eisenstein et al 



and 



Eisenstein et al. 



(120061 ) are reasonable 



fits, b ut in general within 3a of each other, as in 
( 120061 ). We thus conclude that the uncertainties in 
and can now analyze the results without attributing any noted irregularities to the fitting 
pr ocess. 



Kepler et al.l (120061 ). however, compare SDSS spectra with new SNR(g)~ 100 spectra 
acquired with GMOS on the Gemini 8 m telescope for four white dwarf stars around T c g ~ 
12 000 K. Their fits suggest that published SDSS optical spectra fits overestimate the mass 
by AM. ~ 0.13 M.q, because of the correlation between the derived T eff and logg — a 
small increase in T e g can be compensated by a small decrease in logo. Our simulations 
indicate this discrepancy is concentrated only in the region around the Balmer line maximum, 
14000 K^ T cfr ^11000 K. 

To explore the increasing mass trend in more detail, we restricted our sample to the 
1733 stars both brighter than g=19 and hotter than T e g = 12 000 K and obtained an aver- 
age DA mass of (-M)da = 0.593 ± 0.016 M.&. The distribution for this hot and bright 
sample, shown in Fig. [H is s imila r to that of the Palo mar Green survey published by 



Liebert. Bergeron, fc Holbergi (120051 ). They studied a complete sample of 348 DA stars with 
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Table 1. Gaussian fits for the T e ff 12 000 K and g $C 19 histogram, seen in Fig.[T] 



i 


a, 


Mi(M e ) 




Fraction 


1 


264.8 


0.578 


0.047 


69% 


2 


27.8 


0.658 


0.149 


23% 


3 


27.0 


0.381 


0.050 


7% 


1 


3.0 


1.081 


0.143 


1% 



SNR^ 60 spectra and determined atmospheric parameters by spectral fitting via the line 
profile fitting technique, using models up to log g = 9.5. They found a peak in the mass 
histogram at 0.565 Ai & containing 75% of the sample, a low mass peak with 0.403 Ai Q 
containing 10% of the sample, and a high mass peak at 0.780 Ai Q containing 15% of the 
stars. They fit their mass histogram (PG mass histogram from hereafter) with three Gaus- 
sian profiles: 0.565 M Q with a ~ 0.080 M Q , 0.403 M Q with a ~ 0.023 M Q , and a broad 
high-mass component at 0.780 M. & with a ~ 0.108 They found more s t ars aboy e 1 A4 Q 



tha n can be desc r ibed b y the three Gaussians they fit. iVennes et al.l (119971 ). I Vennesl (Il999l ). 
and Marsh et al.l (119971 ) also find an excess of white dwarf stars with masses above 1 M. & 
in their sample of T er r > 23 000 K white dwarf stars. 

The overall mass distribution of our bright sample matches well with that of the previous 
standard PG survey sample. We now explore the distribution with temperature. 



In Fig. [2j we show the mass distribution vs. temperature for DA stars 



brighter than 



g — 19 along with the similar distribution from lLiebert. Bergeron, fc Holbergl (120051 ). Again, 
we see the distributions are roughly equivalent and we see an increase in measured mass at 
lower temperatures. Our histograms use 0.025 M. Q bins (N/dm= const ant) because that is 
the approximately mean uncertainty in our mass determinations. 

To explore the region of increasing mass further, Fig. [3] shows the mass histogram only 
for the 964 DAs brighter than g=19, and 12 000 K ^ T cS ^ 8 500 K, for which we ob- 
tain (M) c ^ 1 ~ 0.789 ± 0.005 M G . We have excluded the stars cooler than T cS = 8500 K 
from our mass histograms because our cooler atmospheric models are not accurate for log g 
determination, as explained earlier. 

Tables [T] and [2] detail the Gaussian fits we made for the histograms of Figures [T] and [3] 
respectively, with 

(M-M t ) 2 ' 



N = J2 a i ex P 



(i) 



Figure 11 of iLiebert. Bergeron, fc Holbergl (120051 ) also shows an increase in mass below 



T er | = 12 000 K (see Fig. [2]), even though they have a limited number of cooler stars due to 
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Sloan DR4 DAs g<19 T trl >12000K <M>=0.593±0.016 M 



'0.57& M s = 69k 
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Figure 1. Histogra m for the 1733 DA stars brigh ter than g=19 and hotter than T e g = 1 2 000 K, compared to the PG surve y 
published bv lLiebert. Bergeron, fc Holberel (|2005h and the SDSS DR1 sample published bv lMadei. Nalezvtv. fc Althaud l|2004l . 
Gaussian fits detailed in Table[T]are also shown. Our bins are 0.025 Mq wide. The second graph shows the DR1 and PG survey 
data normalized to the DR4 sample, even though those samples are smaller and therefore have significantly larger errorbars. 



Table 2. Gaussian fits for the 964 DAs with 12 000 K ^ T eff ^ 8 500 K and g < 19 histogram. 



i 






0i 


Fraction 


1 


29.5 


0.818 


0.160 


49% 


2 


59.6 


0.793 


0.052 


33% 


3 


18.0 


0.640 


0.086 


16% 


4 


13.4 


1.096 


0.136 


2% 




Figure 2. Masses for all 3595 DA white dwarf stars brighter than g=19 and co oler than 40 000 K, showing an i ncreas e in mean 
mass for lower T of j. The large solid (blue) circles are the values published by iLiebert. Bergeron, fc Holbcrg (2005), showing 
the increase in mass at lower T c ff is also present in their sample, which uses a totally independent grid of models and fitting 
technique. 



color se 
used in 



ection effects in the PG survey. It is important to note that the model atmospheres 



Liebert. Bergeron. &: Holbergl (120051 ) and the line profile fitting technique they use 



are totally independent of our own. Therefore, if the observed increase in the measured 
gravity with temperature is merely an artifact of the models, then similar problems must be 
present in two independent groups of models and fitting techniques. We are thus gathering 
increasing evidence that either a) both DAs and DBs really do have higher mean masses at 
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Sloan DR4 DAs g<19 12000K>T eff >8500K <M>=0.789±0.005 

— i 1 1 1 1 1 1 1 1 1 1 1 1 1 — 

0.818 M = 49% 




Figure 3. Histogram for the 964 DA stars brighter than g=19, with 12 000 K ^ T c fj 8 500 K along with the fit Gaussians as 
detailed in Tabled 

lower temperatures, or b) there is a common artifact in the model used for all white dwarf 
spectral fitting. 



4 WHY WE DO NOT TRUST MASSES FOR T EFF < 12000 K FOR DA 
STARS 



Bergeron et al. 



(119951 ) measured an increase in the mean mass for the ZZ Ceti star sample 



around 13000 K> T cS > 11000 K, but indicated it might come from a selection effect because 



the discovery of pulsating stars might have preferred higher mass stars. lArras. Townsley. fc Bildsten 
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(120061 ). e.g., show that the white dwarf pulsators with lower masses should pulsate at cooler 
temperatures. Our sample of 351 bright stars in the same temperature range show a sim- 
ilar increase in mass compared to the hotter sample, but we have not been biased by the 
pulsators, so an observational bias is not the cause for the increase in mass detected. 

The simple expectation that massive stars cool faster than their less massive counter- 
parts does not hold for T c r ^ 10 000 K, as the most massive stars have smaller radius 
and, therefore, their cooling slows down after a few e-folding timescales. Another possible 
explanation for an increased mean mass at lower effective temperatures is the presence of 
otherwise undetected He at the surface, broadening the observ e d H li nes and thus mimicking 
a higher logg. Theoretical models (e.g. iFontaine fc Wesemaell (119911 )) indicate that only for 
DAs with hydrogen layer masses below M.u — 10 _10 .M* will He mix around T e g- = 10 000 K 
and, if A4 h = 10~ 7 .M*, only below ~ 6 500 K. However, all seismologically measure d 



H layer masses are A4h > 10 Ai* (IBradley 



2006, 



2001 



1998 



Fontaine fc Wesemael 



1993). 



Since our increased mass trend happens significantly hotter than T e g- = 6 500 K, He contam- 
ination cannot account for the observed increase in mass at lower temperatures, unless the 



Lawlor and MacDonald's 



more distant stars studied here have significantly thinner H layers 
(120061 ) models show around 3% of the DAs could have M. u ~ 10 _9 .M*, but not thinner. 
Therefore, there are not enough stars with thin H layers, at any rate, to account for our 
exc ess of massive objects. 



Wilsonl (120001 ) proposes a possible physical model for increasing WD masses at lower 



temperatures. She suggests that low metallicity AGBs will produce higher mass white dwarf 
stars, probably around 1M.Q, because the relatively lower mass loss expected for low metal- 
licity AGB stars increases the mass of the core prior to the star moving out of the AGB. 
Since the earlier generations of white dwarf stars which have now cooled more than their 
later cohorts, presumably came from lower-metallicity progenitors, this mechanism could 
explain a mass increase at lower white dwarf temperatures. If we extend this concept to 
globular clusters though, we would expect the mass of the white dwarfs in globular clusters 
to be larger than the mean mass of ou r stars cooler than 10 000 K, which is not observed 



(IMoehler et al. 



2004 



Richer et al. 



20071 ). So again, we are left with a discarded explanation 



of the observed mass increase. 



An interesting clue to the problem may be found in lEngelbrecht fc Koesterl (120071 ) . which 
used SDSS photometry alone to make a mass estimate. Their cool white dwarf stars show 
mean masses similar to those of the hotter stars. Our mass determinations using photometric 
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SDSS DR4 DAs Masses from Colors 
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Figure 4. Masses for all 7167 DA white dwarf stars derived comparing only the SDSS colors (u-g) and (g-r) with those predicted 
from the atmospheric models. For T cf j 20 000K and T g ^ 9 000 K, the colors are degenerate in mass. 



colors only, shown in Fig. HI is derived comparing only the SDSS colors (u-g) and (g-r) with 
those predicted from the atmospheric models convolved with SDSS filters. They do not show 
any increase in mass with decreasing T c g. Because of their larger uncertainties than the 
spectra fitting, we binned the results in 2000 K bins. This result suggests that any problem 
in the models is mainly restricted to the line profiles, not the continua, which dominate the 
broadband photometric colors. 

Thus, we are mainly left with the possibility raised by 



Koester 



( 119911 ) that an increase 



in mass with lower temperatures could be due to the treatment of neutral particles in model 




Figure 5. Masses for all 5718 DAs (crosses) cooler than 40 000 K and 507 DBs (filled circles) showing the continuous increase 
in average mass at lower T c fj ■. 



Bergeron et al. 



(119951 ). however, suggests 



atmospheres with the Hummer-Mihalas formalism, 
that the neutral particles are only important below T e ff ~ 8000 K which is certainly lower 
than where we see the trend begin. It seems the only remaining explanation is that accurate 



modeling of neutral particles will indeed show an effect for DAs near 12 000 K. 



5 DB WHITE DWARFS 

We determined masses for the lEisenstein et al.l (120061) DBs from their fit temperatures and 
gravities using evolutionary grids of lAlthaus et al.l (120051 ); lAlthaus. Serenelli. fc Benvenuto 
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Table 3. Gaussian fits for the histogram of the 150 DBs brighter than g=19 and hotter than T e fj = 16 000 K. 



i 


at 


Mi(M e ) 




Fraction 


1 


8.3 


0.700 


0.109 


59% 


2 


14.3 


0.712 


0.042 


40% 


3 


0.6 


1.288 


0.035 


1% 



( 2001 ). The Althaus mod e 



(120051 ) and 



Metcalfe et al. 



s use time resolved diffusion throughout evolution. 



Metcalfe 



(j2005al ) discuss asteroseismological results in DBs, showing the 
observations are consistent with the layer masses predicted by current diffusion theory. Fig. [5] 
shows we find an increase in the measured surface gravity below T c q ~ 12 000 K for DAs 
and a similar increase below T e g ~ 16 000 K for DBs. For the 208 DBs brighter than g=19, 
we find (M)$ B = 0.785 ± O.O13A^ For the 150 DBs brighter than g=19 and hotter than 
T e s = 16 000 K, we find (-M)db = 0.711 ± O.OO9.M0. Both measurements are considerably 
larger than the 0.593 M. Q mean mass value for the bright and hot DA sample. A similar 
larger (relativ e to th at of the DAs) DB mean mass value has been previously reported by 



Koester et al 



(1200 ll ) who obtained a (-M)db = 0.77 for the 18 DBs they observed with 



UVES/VLT, including stars down to T e g ~ 16 000 K. Others, however, fi n d low er mean 



Oke. Weidemann. fc Koesterl (11984 ) derived 



DB masses, more similar to those of the DAs. 
(M)rm = 0.55±0 .10 from their sample of 25 DBs ranging 30 000 K > T cff > 12 000 K, while 



Beauchampl (I1995T ) found ( ML = 0.59±0.01.M^ for his 46 DBs, ranging 12 000 K ^ T cS ^ 



31 000 K. For the 34 DBs in 



Castanheira et al. 



2006|), ranging 27000 IO T cfr ^ 13 000 K, 



the mean is {M) DB = 0.544 ± O.O57W . 

The Gaussian fits for the 150 DBs brighter than g=19 and hotter than T e g = 16 000 K 
are listed in Table [31 The mass histogram for DBs is shown in Fig. [6j 



6 OBSERVING VOLUME CORRECTION 



In order to turn our observed mass distributions into a real analytical tool, we must first 
correct the sample for completeness. We do this by the l/V max formalism. V max is the volume 
defined by the maximum dista n ce at which a given object would still appear in a magnitude 



limited sample ((Schmidt 



19681 ) 



Geijo et al 



(120061 ) discuss white dwarf luminosity function 
completeness corrections and conclude that for large samples, the l/V max method provides 
a r eliable characterization of the white dwarf luminosity function. 



Liebert. Bergeron. &: Holbergi (120051 ) find that 2.7% of the stars in the PG sample have 



masses larger than 1A4 Q and, when corrected by l/V mauX , 22% are above O.8.M0. 
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Sloan DR4 DB Mass Distribution g<19 T eff >16000K <M> = 0. 71 1 ±0.009 
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Figure 6. Histogram for the 150 DBs brighter than g=19 and hotter than T cS = 16 000 K in DR4. 



We first calculated each star's absolute magnitude from the T e ff anc 



from our fits (for the extreme mass ones) or those of 



Eisenstein et al. 



og q v alues obtained 



(j200J) (for the rest), 



convolving th e synthetic sp ectr a with the g filt e r tran smission curve. We used the evolution- 



ary models of IWoodl (Il995l ) and lAlthaus et al.l (120051 ) with C/O cores up to logg = 9.0, and 
O/Ne cores for higher gravities, Mn e = 10 -2 M*, and Mh = 10~ 4 or M*, to estimate stellar 
radii for DAs and DBs respectively. We do not claim that the SDSS spectroscopic sample is 
complete, but we do contend tha t in te rms of mass, there should be no preferential bias in 

(120061 ) report that spectra are obtained for essentially all 



the target selection. 



Harris et al. 
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white dwarf stars hotter than 22 000 K. Additional white dwarf stars down to T e g = 8000 K 
are also found, but few cool er than that as these sta rs overlap in color space with the F, G, 



and K main- sequence stars. 



Eisenstein et al. 



(120061 ) discuss the spectroscopic sample com- 



pleteness, which is around 60% at 18 < g < 19.5 for stars hotter than T c g = 12 000 K and 
around 40% for cooler stars. Our analysis is restricted to the sample brighter than g=19. 

Once we had our calculated absolute magnitudes, we could estimate each star's distance 
as shown in Fig. 0, neglecting any effects of interstellar extinction. The mean distance for our 
DA samples are: 474 ± 5 pc for the entire 7167 DA sample, 302 ± 5 pc for the stars brighter 
than g=19, and 436 ± 7 pc for the stars brighter than g=19 and hotter than T e fj ~ 12 000 K. 



For each star in our sample, we calculate 

4:71 

V max = y/3(r max - r^ in ) exp(-z/z ) 

where (3 is the fraction of the sky covered, 0.1159 for the DR4 sample, r min is due to the 
bri ght magnitude limit, g=15, an d zp is the di sk scale height whic h we assume to be 250 pc, 



as 



Liebert. Bergeron, fc Holbergi (120051) and 



distribution indicates zq — 310 pc. 



Harris et al 



Vennes et al 



20061 ). even though our height 



(120051 ) show that both the white dwarf 



stars in the SDSS DR1, and the 1934 DAs found in the 2dF (18.25 ^ bj ^ 20.85) quasar 
surveys, belong to the thin disk of our G alaxy. Using these data, they measured a scale 



height around 300 pc. 



Harris et al. 



(120061 ) calculate the white dwarf luminosity function 
from photometric measurements of the white dwarf stars discovered in the SDSS survey up 
to DR3. They assume log g = 8.0 for all stars and use the change in number per magnitude 
bin to calculate the scale height of the disk, obtaining 340l7o° pc, but adopt 250 pc for better 
compa r ison with other studies. This volume inc ludes the disk scale height as discu s sed b y 



Greenl (119801 ) ; iFleming. Liebert. fc Green! (119861 ) ; and iLiebert. Bergeron, fc Holbergi (120051 ) . 



Each star contributes l/V mSuX to the local space density. 

Figures El and [9] show the resulting corrected mass distribution for our DA and DB sample, 
respectively. Figure [8] contains 1733 bright, non-cool DAs, i. e., those with T eS ^ 12 000 K 
and g ^ 19. We also list the corresponding five Gaussian fit parameters in Table HI Figure [9] 
contains 150 bright, non-cool DBs, i. e., those with T c g ^ 16 000 K and g ^ 19. The 
corresponding three Gaussian fits are listed in Table [51 

Since the most massive white dwarf stars have smaller luminosities because of their 
smaller radii, after applying the 1/V max correction to the observed volume, we find that 
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Figure 7. Distribution of distances, d, and height above the Galactic plane, z, for DAs in the SDSS DR4. 



around 20% DAs are more massive than 0.8 Mb in our bright an d hot sample, of the same 
order as that discovered by iLiebert. Bergeron, fc Holbergi (120051 ) for the PG sample. The 
DB distribution is interesting, however, as it tends to significantly higher masses than does 
the DA distribution! We found only two stars from our sample with published atmospheric 
parameters in the litera ture, with AT fi ff = 510 ± 30 and A log q = . 12 ± .15, so we could 



not do a comparison as 
results. 



Kleinman et al 



(120041 ) and 



Eisenstein et al 



(120061 ) did for the DA 
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SDSS DR4 DAs g<19 T eff >12000K 
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Figure 8. Histogram for the 1733 DA stars brighter than g=19 and hotter than T c fj = 12 000 K corrected by 1/V m 



7 EXTREME MASS STARS 



Nalezyty fc Madejl (12004 ) published a catalog of all massive white dwarf stars then known, 
with 112 stars stars more massive than O.8.M0. The four stars with M. ^ 1.3A4 Q in their 
list are ma g netic ones and therefore have large uncertainties in their estimated masses. 

(120041 ) found one non-magnetic massive white dwarf, LHS 4033, with M. ~ 



Dahn et al 



1.318 — 1.335.Mo, depending on the core composition. Our oxygen-neon core mass for their 
derived T eff = 10 900 K and log g = 9.46 is M ^ lSOM^. We not e that the models from the 



Montreal group used to derive T eS and log g in 



Dahn et al. 



(120041 ) show the same increase in 
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Table 4. Gaussian fits for the T e ff 12 000 K and g ^ 19 volume corrected DA mass histogram. 



i 




Mi(M Q ) 


0i 


Fraction 


1 


5.965 x 10" 6 


0.603 


0.081 


38% 


2 


1.203 x 10~ 5 


0.571 


0.034 


32% 


3 


1.165 x 10" 6 


0.775 


0.201 


19% 


i 


1.455 x 10" 6 


1.175 


0.076 


9% 


5 


8.305 x 10" 7 


0.358 


0.037 


2% 



Table 5. Gaussian fits for the volume corrected histogram of the 150 DBs brighter than g=19 and hotter than T e g = 16 000 K. 



i 


a. 


Mi(Me) 




Fraction 


1 


8.6 x 10~ 7 


0.718 


0.111 


66% 


2 


8.9 x 10~ 7 


0.715 


0.045 


27% 


3 


2.9 x 10~ 7 


1.286 


0.031 


7% 



mass with decreasing T e fj as our models and therefore we do not take this mass determination 



as com pletely reliab 
0058), bobbie et all 



'. e due t o the objects relatively low temperature. For GD 50 (WD J0348- 
(12000 ) found T eff = 41550 ± 720 K and logg = 9.15 ± 0.05. Our 
oxygen-neon core mass for their derived T e g and logg is M. — 1.23 ± O.O2.M0, very similar 
to the value reported by them for C/O models. They also show this massive star is consistent 
with its formation and evolution as a single star, not the product of a merger. 

From the 7167 pure DA white dwarf stars, we found 1611 (22%) with A4 > O.8.M0. 
For the 2945 stars brighter than g=19 we found 760 (26%) with M > 0.8.M©, but for 
the 1733 stars brighter than g=19 and hotter than T e g = 12000 K, we find only 105 stars 
(6%) with A4 > O.8.M0. The most massive star in our hot and bright sample is SDSS 
J075916.53+433518.9, whose spectrum spSpec-51883-0436-045 is shown in Fig. QUI with 
g=18.73, T eff = 22100 ± 450 K, log g = 9.62 ± 0.07, M = 1.33 ± O.OIA*©, and estimated 
distance of d = 104 ± 4 pc. We caution that the evolutionary models used to estimate the 
radii, and therefore the masses, in our anal ysis do not include post-newtoni an corrections, 
important for masses above M. — I.3&M0 (IChandrasekhar fc Tooperlll964l ). For the stars 
brighter than g=19, we find 21 others with masses larger than Ai = 1.3 -M , all below T c g = 
9 000 K. We deem the mass determinations for stars cooler than T e ^ ~ 12 000 K unreliable. 
In Table El we list the DAs with 1.2 M & < M < 1.3 A4 hotter than T eff = 12 000 K. 

The spectrum for the brighter g=17.99 SDSS J094655. 94+600623. 4 is shown in Figure fTTl 
Because our analysis uses relatively low SNR spectra and gravity effects dominate mainly 
below 3800A, where we have no data, our conclusion is that we must undertake a study 
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SDSS DR4 DB Mass Distribution 
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Figure 9. Histogram for the 150 DBs brighter than g=19, and hotter than T c fj = 16 000 K in DR4, corrected by volume. 



Spectra-M-P-F 


Name 


g 


M g 


T e ff 


(7T 


logs 




M 


0"A1 


d 










(K) 


(K) 






(M e ) 


(Me) 


(pc) 


spSpec-51691-0342-639 


SDSS J155238.21+003910.3 


18.44 


12.23 


15924 


387 


9.280 
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1.262 


0.010 


97 


spSpec-51915-0453-540 


SDSS J094655.94+600623.4 


17.99 


10.87 


28125 


220 


9.370 


0.040 


1.287 


0.010 


123 


spSpec-52374-0853-198 


SDSS J133420.97+041751.1 


18.52 


12.34 


17549 


422 


9.150 


0.060 


1.223 


0.020 


125 


spSpec-52703-1165-306 


SDSS J150409.88+513729.1 


18.84 


10.28 


79873 


8228 


9.050 


0.390 


1.204 


0.100 


468 


spSpec-52751-1221-177 


SDSS J110735.32+085924.5 


18.42 


12.23 


18715 


327 


9.140 


0.060 


1.219 


0.020 


128 


spSpec-52872-1402-145 


SDSS J154305.67+343223.6 


18.33 


10.85 


30472 


313 


9.300 


0.070 


1.269 


0.010 


168 



Table 6. DA stars with masses above 1.2 Me ana - t> e l° w 1-3 Me derived from the SDSS spectra, with T c ff 12 000 K. 
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SDSS J075916. 53+433518.9 g=18.73 T ef( =22100±450K M=1.33±0.01 M 




4000 5000 6000 7000 

MA) 

Figure 10. Spectrum of SDSS J075916.53+433518.9 with g=18.73, T off = 22100 K and two models, with logg = 9 and 10. 
The higher log g fits the Ha line better, but the lower log g fits the higher lines better, where the SNR is smaller. This low SNR 
is typical for the stars closer to our upper cutoff of g=19. 

in the violet or ultraviolet to measure the masses more accurately. An extensive study of 
gravitational redshift would also be critical. 

For the 507 single DBs we find 30 DBs with log g > 9. Most of our massive DBs are 
cooler than T eff ~ 16 000 K, or fainter than g=19, except for SDSS J213103.39+105956.1 
with g=18.80, T cS = 16476 ± 382, and logg = 9.64 ± 0.21, corresponding to a mass M. = 
1.33 ± 0.04.M©, and for SDSS J224027. 11-005945.5 with g=18.82, T cS = 17260 ± 402, and 
logg = 9.31 ± 0.20, corresponding to a mass Ai = 1.25 ± 0.067^©. 

The low mass stars present in our sample are consistent with He core evolution models 
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SDSS J094655. 94+600623.4 g=18 T eff =28100±220K M=1.29±0.01 M 
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Figure 11. Spectrum of SDSS J094655. 94+600623.4 with g=17.99, T cS = 28100 K and two models, with logg = 9 and 9.5. 
This spectrum is typical of the the SNR achieved for the 1003 DAs and 59 DBs brighter than g=18 in our sample. 



calculated by 



Althaus. Serenelli. fc Benvenutol (120011 ). and displayed in Fig{T21 It is impor- 



tant to stress that these stars should be studied with more accurate spectra and model 
atmospheres, as they are possible progenitors of SN la if they accrete mass from compan- 
ions. 



8 CONCLUSIONS 

Our investigatio n s into the mass distribution of the SDSS DR4 white dwarf sample from 
Eisenstein et al.l ( 120061 ) revealed several items. First, all groups are seeing nearly identical 
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Figure 12. Evolutionary tracks for He white dwarf stars calculated bv lAlthaus. Serenelli. fc" Benvenutd |200lft and the location 
of the lowest mass stars in our sample. Most of these stars are below g=19 and therefore have noisy spectra. 



increases in mean DA at lower temperatures (less than « 12 000K for DAs and 16 000K for 
DBs). Either this is truly going on in the white dwarf stars, or there is missing or incorrect 
physics in everyone's models. We propose the treatment of neutral particles as the most 
likely explanation. We suspect the atmospheric models should be improved with a detailed 
inclusion of the line broadening by neutral particles, since the increase in apparent mass for 
both DAs and DBs occur at temperatures when recombination becomes important. 

Secondly, we find a significant difference between the DA and DB mass distributions, with 
the DB distribution significantly more weighted to massive stars. Figs. [T3] and show the 
combined number and volume-corrected DA and DB histograms. The DB histograms have 
been re-normalized to the DA maximum for display purposes. Our results contradict nearly 



al 
of 



previous work whic h show the mean DA and DB masses to be similar (with the exception 



Koester et al 



200 ll ). We note that the previous efforts, though, were based on histograms 



for DBs with less than 50 stars, and our DB histogram has 150 stars. However, we still need 
to explore our DB models and fits in more detail to verify the validity of this novel result. 
Specifically, we find (M) DB = 0.711 ± O.OO9A^ , higher than (M) DA = 0.593 ± 0M6M Q 



White Dwarf Mass Distribution in the SDSS 23 
Sloan DR4 



300 



200 - 



100 - 




0.2 0.4 



0.6 0.8 1 

Mass (M ) 



1.2 1.4 



Figure 13. DA and DB histograms for comparison. The DB histogram has been re-normalized to the DA maximum for display 
purposes. 



for the 1733 DAs brighter than g=19, and hotter than T e g = 12000 K. This is a significant 
new result and must be investigated further. 

We have also detected a large number of massive DA white dwarf stars: 760 with M. > 
O.8A^0 brighter than g=19 and 105 both brighter than g=19 and hotter than T e s = 12000 K. 
We report the highest logg white dwarf ever detected. 
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Figure 14. DA and DB histograms corrected by observed volume for comparison. The DB histogram has been re-normalized 
to the DA maximum for display. 



